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Abstract 

In 1+1 dimensional hydrodynamics originally proposed by Landau, we derive a new poten- 
tial and distribution function including Heaviside function and investigate its mathematical and 
physical properties. Using the original distribution derived by Landau, a distribution function 
found by Srivastava et al., our distribution function, and the Gaussian distribution proposed 
by Carruthers et al, we analyze the data of the rapidity distribution on charged pions and K 
mesons at RHIC energies ( ^/sjvat = 62.4 GeV and 200 GeV ). Three distributions derived from 
the hydrodynamics show almost the same chi-squared values provided the CERN MINUIT is 
used. We know that our calculations of hadron's distribution do not strongly depend on the 
range of integration of fluid rapidity, contrary to that of Srivastava et al. Finally the roles of the 
' Heaviside function in concrete analyses of data are investigated. 

in ' 
in 

CO 



00 

o 



PACS. 25.75.-q Relativistic heavy-ion collisions - 24.10.Nz Hydrodynamic models 



1 Introduction 



Recently authors of ref. [TJ [2] have mentioned a paper by Srivastava et al. [3] as one of boost non- 
invariant solutions in 1 + 1 dimensional hydrodynamics proposed by Landau [U [6] . It is well know 
^ ■ that the boost invariant distribution function have been proposed by Hwa [7] and independently 
Bjorken [8]. It is also worthwhile to notice recently interesting papers on 1 + 1 dimensional hydro- 
dynamics published in refs. [9l[T0]. See also ref. [TT] in which the usefulness of the hydrodynamics 
for the collisions (proton-proton and nuclei- nuclei) is stressed. 

In the present study, in particular we are interested in the paper by Srivastava et al., because 
their potential and the distribution function are similar to that calculated by Landau and coauthors. 
Our purpose is to obtain a new solution by use of different boundary conditions in the same partially 
differential equation obtained by Landau. (See eq. (jBJ) below) 

The hydrodynamics is briefly described in the following: The energy density e, the pressure 
p, the four velocity u^, and the entropy density s are necessary physical quantities. Then the 
conservation law of the energy-momentum and that of the entropy are given as 
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where T^ v = (e + p)u^u v — pg^ u . Introducing the temperature of the fluid T, the equation of state 
p = (? s e with the velocity of sound c s , the thermodynamical relations de = Tds, and e + p = Ts, 
we have the relation dp = sdT. Moreover, introducing the rapidity of fluid y and x^ 1 = (t, x) = 
(t cosh r/, rsinhr/) with r\ = tanh _1 (x/t), 1 + 1 dimensional hydrodynamics in the perfect fluid is 
described as follows 

— (rsinhy) + — (Tcoshy) = 0, (2) 

— + c s - ) cosh y + + c s - ) smh „ = 0. (3) 

Utilizing the Legendre transformation d\ = d{(j> + tT cosh y — xTsinhy) from the variables (t, x) 
to (T, y) with potentials (f> and Xi we obtain the following expressions: 

d~x_ 1 d\ 

t = -^-z; cosh y — — — — sinh y, (4) 
T dy K ' 

()y 1 5% 

x =-— sinh w — ——— cosh -u. (5) 
dT y Tdy y w 

Moreover, we have the following relation between y and rj: 



y = tanh" 1 (?) + tanh" 1 (^EL \ = 77 + A. (7) 

From the above equations and relations, we obtain the following partial differential equation of the 
potential x in 1 + 1 dimensional hydrodynamics, 

where uj = — ln(T /Tq) (the logarithmic temperature). To look for a new solution in eq. ([8]) is our 
aim, as mentioned before. 

As the solution of eq. ([8]) , we have known the following two solutions [U [3] . (Notice that two 
different initial conditions are used.) Utilizing the method of Green function |12|, I14j . we obtain 
a new potential including the Heaviside function for the differential equation in 1+1 dimensional 
hydrodynamics . 

On the other hand, the rapidity distributions of the charged mesons at RHIC energies (y / ijv r /v = 
62.4 GeV [E] and 200 GeV [T6]) have been reported by BRAHMS Collaboration. We are going to 
analyze them by several formulae by Landau, Srivastava et al, and ours. 

In the second sec. we briefly explain the formula by Landau and the boost non-invariant formula 
calculated by Srivastava et al. In the third sec. we present the new formula including the Heaviside 
function. In the 4th one we present analyses of the data at RHIC energies by three formulae and 
a formula proposed by Carruthers et al. |17j in addition to them. In the final one the concluding 
remarks are presented. 



2 Solutions of Potential (eq. ([8])) by Landau and Srivastava et al. 



We mention the solution by Landau and that by Srivastava et al. An original formula by Landau 
is obtained by the following initial conditions, 

dx 



dy 



Ule u 



x(y = uj/c s ) = 0, 



(9) 
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where I denotes the length of colliding matter. 

The potential is expressed by the integral formula 



X ~ e 



where (3 = (1 — c^)/2c^. The entropy distribution is calculated as follows, 

dS dN d ( d X 

-r- ~ const — x + -3- 

ay ay acj \ ouj 

where Uf = — ln(Ty/To). Using (fTTj) . we have the rapidity distribution 



dN 
dy 



~ e 



/ (p) + ^/i(p) 



(10) 



(11) 



(12) 



where p = (3\Juj 2 — c 2 y 2 . Iq and I\ are the modified Bessel functions of the 1st order and 2nd one, 
respectively. 

On the other hand, Srivastava et al. have assumed that X = Xi exp(/3u;) and obtained the 
following equation 



9 2 Xl 2^X1 _ 2n2 

dy 2 ° s doj 2 ° sP Xl ' 
Using a new variable with p = (3\/uj 2 — cjy 2 , 

-1 ( c s y 



q = tanh 



UJ 



(13) 



(14) 



they obtain the following equation which is satisfied with the modified Bessel function in a special 
case (as below), 



1 d 2 xi d 2 xi Idxi 



p2 Qg2 Qp2 p Qp 



+ xi = o. 



(15) 



Adopting the condition = and assuming the finiteness on the characteristic (u) = ±c s y), they 
have obtained the following solution, 



X 



e^IoiP^u 2 - cjy 2 ) for uj 2 > cjy 2 , 
e^JoiP^/c^y 2 - to 2 ) for J 2 < cjy 2 . 



(16) 



To obtain the potential with the Bessel function, they have introduced a different variable, q' = 
tanh" 1 (uj I c s y) , in their calculation. To show some mathematical properties of two solutions, we 
present their contour maps in fig. [TJ 

The entropy distribution related to eq. (fl~6j) is calculated as 



dS 
dy 



P 



(17) 



for uj > cjy . For uj < cjy we have 



dS_ 

dy 



p* 



(18) 



where p = Py/c'jy 2 - uj 2 . Authors of refs. [H EE [II] have applied eqs. (UJ), (HZD and (TTH]) to 
analyses of data. 
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Landau, cJ=l/3 
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(e) yc, 

Figure 1: Distributions and contour maps of % °f ec l s - POP an d (I16p . (a) Distribution of eq. (jlQI) 
with = 1/3. (b) Contour map of fig. Ufa). ( c ) Distribution of eq. (fl~6|) with <~ = 1/3. (d) Contour 
map of fig. ttfc). (e) When assumption with x( w > 2/) = for w < c s |y| is used, this figure is obtained. 
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3 Solution including the Heaviside function 



Hereafter, to consider a new solution, we use the method of Green functions [12\ [T3"l Q3] 



d 2 X 
duj 2 



dx 1 9 2 x 



dco c 2 dy 2 



-^5(uj - uj )5(y - y ). (u , y -> 0) 



(19) 



Using the function x = Xi e > we choose the following initial conditions in the (u>, y) space, because 
we can consider the initial conditions at u = 0, 



dxi 



LO = 0+ 



G{y) = 6(y), X i(w = 0+, y) = g(y) = 0. 



(20) 



The reason is as follows: In the localized point ((o>, y) ~ (0, 0)), the potential should exist as 
the boundary value (as the initial condition), then we choose 5(y), because the fluid created in 
collisions is localized there. Moreover, the initial condition of the potential at u = is chosen as 
x{w = 0+, y) = 0. See eq. ([9]). This set of initial conditions suggests us that the present approach 
mainly describes the collision in the central region. 

In order to obtain the potential governing the hydro-system in collisions, we perform the following 
calculation named the Riemann's formula, 



Xi(y,v) 



^{g(y + u/c a ) + g(y - uj/c s )} 



+ 



~ \c s G(z + y)h{pJuj 2 - c 2 z 2 ) + (3 2 c s u;g(z + y) 

£ J-u)/c a y 



h((3^ 2 -cfz 2 ) 

flyjuj 2 — cjz 2 



Using eq. (I20p . we have the following solution 

X ~ e^I (pJu; 2 -c 2 y 2 ) H(Q = lo - c s \y\), 



>dz. (21) 



(22) 

where H{Q) is the Heaviside function defined as follows 

ff<g-»-*M)-{j q<S: 

The Heaviside function in eq. (|22|) appears from the region of integration in eq. (|21[) , Hereafter, in 
concrete analyses, introducing a parameter e, we use the following approximate expression, 



H(e, Q) 



1 



+ 1 



(23) 



To show some mathematical properties of the above solution, we present its contour map in fig. [2l 
Moreover, we show diagrams between (t, x) for eq. (fT6|) and ([22]) in fig. [3] (a), (b). In fig. [3] (a) we 
confirm calculations in ref. [3], and in fig. [3] (b) our calculation for eq. ([22]) is presented. 
Here we examine an evolution of the temperature as, 



T 



1 


' ( dx) 2 _ ( dx) 2 ' 




e 2uj 


'(d x 


) 2 -(%) 2 ; 




l 


\dw J \oy J 


u=0, y=0 \ 




m 


2 -( 


dy) 


u=Q, y=0 



(24) 



and 



T 3 to 



—3u> 



TQ 



(25) 



5 



c?=l/3, e=5 
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(c) y°s 

Figure 2: (a) (b) Distribution and contour map of eq. (I22D with e = 5. Behaviors distinguished 
from that of fig. [I] are observed along the lines u = ±yc s . (c) Contour map of eq. ([22]) with e = 50. 
The figure is the same as fig. [He) as e = oo. 
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Figure 3: Trajectries of t = \/r 2 + x 2 (r = 1 fm). (a) Boost invariant solution by Hwa and Bjorken 
(dashed lines) and boost non-invariant one by Srivastava et al. (solid line), (b) Boost invariant 
solution (dashed lines) and our boost non-invariant one by us with e = 2.5 (solid line). Finite 
difference method is used. 
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Figure 4: The temperatures space-time evolution at y = and y = 2 by means of eq. (j25[) are 
displayed. Notice that results at y = by Srivastava e£ aZ. and ours with small e = 2.5 are the 
almost same. Ours with e = 5 is similar to result by Landau. At y = 2, we cannot estimate a figure 
by Landau's original solution, because of constraints in integrations. 
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Figure 5: (a) Behavior of eq. (f23l) . (b) Behavior of eq. (|26|) . 



The results at y = and y = 2 are presented in fig. 01 We can compare differences of the tempera- 
tures evolution among refs. [31 HI [H] and ours. 

The following expression is the entropy distribution, 



dN 
dy 



dS_ 

dy 



P(J3 + l)I (p) + 



p 2 ( P 2 +p 2 (2(3 + l)u - P 2 u 2 )r ( P ) 0WI%( P ) 



+ 



(2(3 + l)J (p) + 



2f3 2 tol (p) 



P 



H'(s, Q)+I (p)H"(e, Q). 



H(e, Q) 



(26) 



When p is imaginary, the modified Bessel function should be replaced by the Bessel functions, 
because we treat e as a free parameter in the CERN MINUIT program. This means that we are 
able to obtain any useful information for collisions, provided that there was a relation among e, uj 
and c s governed collisions. The typical behavior of eq. (|26p is shown in fig. [5j 

To analyze the measured distribution of hadrons, the following calculation is necessary, 



dN 
dyh 



yo 

dy 
-yo Jo 



dprPT 



m T cosh(y fe - y) 



dN 



exp[{m T cosh(y h - y) - hb}/T] + 5 dy ' 



(27) 



where yh, Pt^t = yPr + m p an< ^ ar e the rapidity of hadron, the transverse momentum of 
hadron and the baryonic chemical potential of hadron A, respectively. For mesons, 5 = —1, [ib = 
and for baryons 5 = 1, = ^- As you see fig. El due to the derivative of the Heaviside function, 
there are small-finite contributions outside of the asymptotic lines uj = ±c s |y|. 



4 Analyses of data by several formulae 



Since BRAHMS Collaboration has reported the rapidity distributions of charged n mesons and K 
mesons at RHIC -^/saF/v = 62.4 GeV and 200 GeV, we can analyze them by three formulae mentioned 
above. The freeze-out temperature and chemical potential are cited from ref. [21]. Tf = 0.1605 (200 
GeV and 62.4GeV) and // = 0. See also ref. [22]. Moreover, we can analyze the data by the Gaussian 
distribution proposed by Carruthers and Duong- Van [17J. 



dN 
dyh 



2tto~ 



cxp 



2a 2 



In 



2m n 



(28) 
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Our results by eq. (|27H in which c 2 s is assumed as a free parameter are presented in figs. [6] and [71 
Notice that we examine different cutoff values for uq in calculations of eq. (|27p . Analyses of data 
by eq. (|28|) are cited in tableland shown in figs. and [7J 



Table 1: Analysis of charged pions and K mesons rapidity distributions by means of eqs. (|12j) and 
(|27j) with yo = u/i/c$. (Distribution by Landau) 





data 


UJ 


c 2 


c 


X7n.d.f 


7T 


(200 GeV) 


1.77±0.03 


0.194±0.016 


1230±426 


18/11 


7T+ 


(200 GeV) 


1.84±0.02 


0.135±0.001 


169±12 


10/11 


K~ 


(200 GeV) 


1.58±0.05 


0.179±0.024 


691±371 


6.3/9 




' (200 GeV) 


2.49±0.02 


0.066±0.001 


0.0013±0.0002 


1.8/9 




(62.4 GeV) 


1.61±0.01 


0.333±0.059 


4860±92 


13/7 


7T+ 


(62.4 GeV) 


1.49±0.16 


0.295±0.074 


4250±1580 


15/7 


K~ 


(62.4 GeV) 


1.22±0.23 


0.189±0.107 


1040±1250 


1.3/5 


K+ 


(62.4 GeV) 


1.40±0.08 


0.178±0.043 


724±503 


0.44/5 



Table 2: Analysis of charged pions and K mesons rapidity distributions by means of eqs. (|17|) . (1180 
and (|27p with yo = w /\/c? + 8 (5 = 2). (Distribution by Srivastava et al.) 





data 




c 2 


c 


X 2 /n.d.f 


7T 


(200 GeV) 


1.09±0.15 


0.246±0.024 


1140±659 


17/11 


7T+ 


(200 GeV) 


1.73±0.02 


0.128±0.001 


6.52±0.46 


10/11 


K~ 


(200 GeV) 


0.90±0.20 


0.236±0.032 


657±500 


6.4/9 


K+ (200 GeV) 


2.24±0.04 


0.102±0.001 


0.0524±0.0090 


1.8/9 


7T 


(62.4 GeV) 


0.48±0.02 


0.333±0.033 


4120±77 


14/7 


7T+ 


(62.4 GeV) 


0.43±0.02 


0.333±0.245 


4200±211 


15/7 


K~ 


(62.4 GeV) 


0.41±0.13 


0.244±0.098 


839±1100 


1.3/5 


K+ 


(62.4 GeV) 


0.66±0.12 


0.240±0.039 


703±464 


0.42/5 



5 Concluding remarks 

We summarize the present studies as follows: 

1) We have investigated the new solution including the Heaviside function in 1+1 dimensional 
hydrodynamics. The Heaviside function is approximately expressed by eq. (|23p . Its typical behaviors 
are shown in fig. [SJ The concrete potential and the entropy distribution are presented in fig. [SJ 
Moreover, three concrete components of our distribution are shown in fig. [H Of course, the second 
and third ones depend on the values of e's. The role of the parameter seems to be reflecting the 
temperature fluctuation on the characteristic (u> = ±c s y) in the present analyses. 

2) Three distributions by Landau, Srivastava et al, and ours show almost the same minimum chi- 
squared as CERN MINUIT is used. However, values of u's are different, provided that eq. ([26]) are 
applied to the data. Eq. ([27]) by Srivastava et al. strongly depends on cutoffs values yo's at 62.4 
GeV. In other words, for values of cutoff no principle exists. The main reason is attributed to their 
assumption of the finiteness on the characteristic (u = ±c s y). Please compare fig. [TJe) and fig.EJc). 

3) On the other hand, our distribution function, due to the Heaviside function, does not strongly 
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Table 3: Analysis of charged pions and K mesons rapidity distributions by means of eqs. (I17h . (1181) 
and (|27p with yo = ln(y / iiV7v/m p ). (Distribution by Srivastava et al.) 









r 2 

c s 


c 


v 2 In c\ f 


7T 


(200 GeV) 


87±0 68 


268±0 065 


2030±2910 


17/11 


7T+ 


(200 GeV) 


1.73±0.02 


0.128±0.001 


6.51±0.46 


10/11 


K~ 


(200 GeV) 


0.69±1.59 


0.255±0.259 


1100±5110 


6.5/9 




(200 GeV) 


2.24±0.04 


0.102±0.001 


0.0542±0.0093 


1.8/9 




(62.4 GeV) 


4.12xl0~ 7 ±0.35 


0.265±0.002 


3470±62 


21/7 


7T+ 


(62.4 GeV) 


1.40xl0~ 6 ±0.60 


0.259±0.002 


3230±59 


18/7 




(62.4 GeV) 


4.55xl0~ 5 ±9.84 


0.247±0.011 


1330±209 


1.4/5 


K+ 


(62.4 GeV) 


1.17xl0~ 6 ±6.47 


0.283±0.008 


2120±134 


0.56/5 



Table 4: Analysis of charged pions and K mesons rapidity distributions by means of eqs. (I26p and 
(1271) with y = uj/y^ + 5 (5 = 2). (Our distribution) 





data 


UJ 




c 


e 


X 2 /n.d.f 


7T 


(200 GeV) 


1.55±0.15 


0.276±0.052 


1120±1190 


1.71±0.44 


16/10 


7T+ 


(200 GeV) 


1.86±0.02 


0.112±0.001 


1.35±0.09 


3.35±0.80 


10/10 


K~ 


(200 GeV) 


1.40±0.14 


0.165±0.010 


50.9±18.6 


6.15±7.66 


6.6/8 


K+ 


(200 GeV) 


2.28±0.04 


0.105±0.001 


0.064±0.010 


2.50±1.61 


1.9/8 




(62.4 GeV) 


0.86±0.02 


0.333±0.331 


2870±80 


3.22±0.43 


8.5/6 


7T+ 


(62.4 GeV) 


0.77±0.03 


0.333±0.223 


3100±140 


2.52±0.33 


16/6 


K~ 


(62.4 GeV) 


0.68±0.09 


0.333±0.267 


1570±220 


2.12±0.74 


1.2/4 


K+ 


(62.4 GeV) 


1.16±0.18 


0.217±0.132 


241±822 


3.66±3.86 


0.40/4 



Table 5: Analysis of charged pions and K mesons rapidity distributions by means of eqs. (I26p and 
(f2*7|) with uq = ln(y / J/v]v/mp). (Our distribution) 



data 


UJ 


c 2 


c 


e 


X 2 /n-d.f 


7T - (200 GeV) 


1.55±0.09 


0.277±0.030 


1150±690 


1.70±0.32 


16/10 


vr+ (200 GeV) 


1.86±0.02 


0.113±0.001 


1.35±0.09 


3.35±0.80 


10/10 


K~ (200 GeV) 


1.37±0.16 


0.167±0.007 


58.9±28.3 


7.31±6.82 


6.5/8 


K+ (200 GeV) 


2.28±0.04 


0.105±0.001 


0.0639±0.0098 


2.50±1.61 


1.9/8 


7T _ (62.4 GeV) 


0.90±0.08 


0.306±0.047 


2150±1210 


3.86±1.35 


8.6/6 


tt + (62.4 GeV) 


0.93±0.02 


0.227±0.026 


741±334 


6.31±3.12 


14/6 


K- (62.4 GeV) 


0.89±0.10 


0.215±0.076 


291±439 


4.95±3.73 


1.1/4 


K+ (62.4 GeV) 


1.16±0.18 


0.217±0.246 


243±857 


3.65±3.93 


0.40/4 
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Table 6: Analysis of charged pions and K mesons rapidity distributions by means of eq. (|28|) . 
(Distribution by Carruthers et al.) 



data 


(7 


c 


xV n -d-f 


7T- (200 GeV) 


2.16 


1626±8 


101/13 


7T+ (200 GeV) 


2.16 


1621±8 


45/13 


K~~ (200 GeV) 


2.16 


244±3 


9.2/11 


K+ (200 GeV) 


2.16 


264±3 


24/11 


7T- (62.4 GeV) 


1.87 


1013±11 


82/9 


vr+ (62.4 GeV) 


1.87 


986±10 


117/9 


K- (62.4 GeV) 


1.87 


130±4 


14/7 


(62.4 GeV) 


1.87 


161±4 


4.2/7 




Figure 6: Comparison of pions and K mesons rapidity distribution and several formulae at y / ijv r /v 
200 GeV. Results by Srivastava and ours in tables [2] and 2] are used. 
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Figure 7: Comparison of pions and K mesons rapidity distribution and several formulae at y / sjv r /v = 
62.4 GeV. Results by Srivastava and ours in tables [2] and 0] are used. 

depend on cutoffs at both energies. It is interesting that three sound velocities c s 's at 62.4 GeV 
are 0.333, and decrease at 200 GeV in tabled! Among the Gaussian distribution and distributions 
based on 1+1 dimensional hydrodynamics there are significant differences as shown in figs. [6] and [71 

4) In conclusion, it can be added that the temperatures in collisions are estimated as lo (200 GeV) 
> oj (62.4 GeV): T = 0.7 GeV (200 GeV) and T = 0.4 GeV (62.4 GeV), respectively, provided that 
the explicit transition (hadrons 44> QCD) (See [19]) was not taken into account. In a future we are 
going to investigate the role of the phase transition between hadrons and QCD in this scheme [19] 
and problem of the fluctuation |23| . 
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